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Abstract 

We consider the quark-anti-quark potential on the three sphere or the generalized cusp 
anomalous dimension in planar M = 4 SYM. We concentrate on the vacuum potential in the 
near BPS limit with L units of R-charge. Equivalently, we study the anomalous dimension 
of a super- Wilson loop with L local fields inserted at a cusp. The system is described by a 
recently proposed infinite set of non-linear integral equations of the Thermodynamic Bethe 
Ansatz (TBA) type. That system of TBA equations is very similar to the one of the spectral 
problem but simplify a bit in the near BPS limit. Using techniques based on the Y-system 
of functional equations we first reduced the infinite system of TBA equations to a Finite set 
of Nonlinear Integral Equations (FiNLIE). Then we solve the FiNLIE system analytically, 
obtaining a simple analytic result for the potential! Surprisingly, we find that the system 
has equivalent descriptions in terms of an effective Baxter equation and in terms of a matrix 
model. At L — 0, our result matches the one obtained before using Localization techniques. 
At all other L's, the result in new. Having a new parameter, L, allows us to take the large L 
classical limit. We use the Matrix Model description to solve the classical limit and match 
the result with a string theory computation. Moreover, we find that the classical string 
algebraic curve matches the algebraic curve arising from the matrix model. 
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1 Introduction 

Recent developments in quantum field theory suggest that for the first time we may be able to 
solve an interacting gauge theory in four dimensions. Most of these developments have come 
about in the study of M = 4 SYM in its planar limit. Planar observables in that theory are 
mapped to two dimensional problems that are solved exactly using Integrability methods 1 . 
Such observables include the spectrum, Wilson loops, scattering amplitudes and correlation 
functions. Among these, the most developed study so far is of the quantum spectrum of the 
theory [4]. It is now available in a completely non-perturbative fashion numerically or even 
analytically in some cases [5-10]. 

Nevertheless, the solution of the spectral problem remains to be rather involved tech- 
nically [6-9] and even sometimes scary [11, 12]. It is given by an infinite set of functional 
equations (Y-system) supplemented with definite analytical properties [13,14]. 

Recently, the methods developed for the spectrum of the theory on S 3 were extended 
to a new class of observables - the spectrum of the color flux between external quarks on 
S 3 [15, 16]. In its vacuum, the energy of that flux is the same as the generalized cusp 
anomalous dimension r cusp . That is, the conformal dimension of a quark and anti-quark 
lines meeting at a cusp [17] 



Here W is a cusped Wilson loop representing the quark-anti-quark lines and Auv, Air are 
the UV and IR cutoffs. 2 The solution takes the form of a Thermodynamic Bethe Ansatz 

1 The integrability in QCD was first discovered in [1,2]. Very fast development of the integrability in 
N = 4 SYM started after the seminal paper [3] . 

2 The phrase "cusp anomalous dimension" is also used to describe the IR divergence of massless particles. 
That cusp anomalous dimension T^ sp is one of the first observables computed exactly using integrability [5] . 
It controls the behavior of T cusp in the limit where the two quarks are infinitely boosted with respect to each 




(1) 
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(TBA) system very similar to the one for the spectrum. These are an infinite set of non-linear 
integral equations. 




Figure 1: (a) A quark- -anti- quark pair sitting at two points on S 3 at a relative angle tt — <fi. 
The quark-anti- quark lines are extended along the (Euclidian) time direction, (b) Under the 
cylinder to plane conformal map, the quark and anti-quark lines in (a) are mapped to the 
two half lines of a Wilson line with a cusp angle (ft. 



The most natural quarks in = 4 SYM are infinitely massive w-bosons on the boundary 
of the Coulomb branch. These are locally supersymmetric quarks probes that also couple to 
a scalar. As we have six scalars, this coupling selects a point n in S 5 . As a result, one of 
the new key features of the cusp TBA system is that it is parameterized by two continuous 
parameters. That is, r cusp is a function of two angles (ft and 9 [18]. The angle is the 
geometrical angel between the two lines, see figure 1. The second angel 9, is the angel on S 5 
between the quark and anti-quark points cos 9 = n q ■ tig. 

The existence of these two parameters allows one to take various limits where the TBA 
equations simplify. We hope that a better handle on the solution in such limits will shed 
light on how to simplify the TBA equations in general. One such limit is (f 2 — 9 2 — > 0. When 
<f) 2 = 9 2 the cusped Wilson loop is BPS and the energy vanishes [19]. Consider the case 
9 = 0. As we deform the angles away from this supersymmetric configuration, the energy 
behaves as 

T cusp (<j ) ,Q) = - ( j ) 2 B (g) + O(<f ) A ) . (2) 

The function Bq, also known as the "Bremsstrahlung function", is related to a variety of 
physical quantities [20,21]. In particular, it controls the power radiation of an accelerating 
quark. That function was computed exactly in [20,21] using results from localization [26,27]. 

In the planar limit the localization result gives 
other lim r cusp = i<f>/4 x We hope that this will not cause confusion. 

i<ft— too ' 
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where I n are modified Bessel functions. On the other hand, the same function B is also 
computed by the TBA. In the near BPS limit the TBA reduces to a somewhat simplified 
set of equations [15], as we review in the next section. These are still an infinite set of 
non-linear integral equations. The simplicity of (3) strongly suggests that, in the near BPS 
limit, the TBA equations can be drastically simplified. Indeed, we managed to solve the 
TBA equations analytically in this limit and derive (3). 

Furthermore, one can consider a more general configuration with adjoint fields inserted at 
the cusp. 3 The addition of L complex scalars Z at a BPS cusp preserves supersymmetry [22]. 4 
That freedom allow us to generalize the Bremsstrahlung function Bq to the case with Z L 
insertion - Bl. We solve the TBA for Bl analytically at any L\ The result is again a rational 
function of modified Bessel functions I n (4:irg): 
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The paper is organized as follows. In section 3 we first apply the methods of [14, 23] 6 to 
reduce this infinite set of equations to just three simple integral equations. Then, in section 4, 
we solve these three equations, obtaining an exact analytic result at any value of the 't Hooft 
coupling, therefore establishing the drastic simplification of the TBA system! The result of 
course agrees precisely with (3) [20,21]. In section 4 we analytically solve the TBA system 
with L units of R-charge, obtaining our main result (4). Very likely, the same localization 
techniques used to derive (3) [20,21,26,27], can be extended to this case. Therefore, our 
results provide for the first time predictions for a localization computation from integrability. 
In section 5 we consider the classical large L limit of the construction. To take this limit we 
were forced to reformulate our result in terms of a matrix model partition function which 
we solve in the classical (thermodynamic) limit. The result is shown to match the classical 
string prediction derived in appendix E. 



2 Bremsstrahlung TBA 

In this section we review the cusp TBA system and its reduction in the near BPS limit 
02 _ q2 _^ name iy the Bremsstrahlung TBA [15]. We also use it to introduce our notations. 

3 In particular, the TBA equations where derived by first considering the problem with infinite fields 
insertions at the cusp. 

4 Here Z is a complex scalar that is not made of any of the scalars coupled to the Wilson lines. 

5 The dashed line shows where a row was deleted. 

6 For a more recent alternative approach see also [24,25]. 
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For simplicity, we will restrict our discussion to the case where 9 = and therefore is the 
smallest parameter.' A reader familiar with [15] may jump to the next section. 

The TBA equations for the spectrum of the color flux between a quark and an anti- 
quark are very similar to the ones describing the spectrum of single trace operators. I.e. 
they follow the pattern of the T-hook with three wings. More precisely, the TBA equations 
divided by the asymptotic solution take exactly the same form. The only differences are in 
the asymptotic solution and a projection 8 . That projection is irrelevant for the vacuum in 
which all Y-functions are even. 

One important feature of the momentum carrying Y-functions is that they have a double 
pole at zero mirror momentum for any value of the coupling, associated with the exchange 
of a single particle [28]. In the near BPS limit all the dynamics of the momentum carrying 
Y-functions reduces to the residue of that pole C a . That is 

\ua(u 2 Y afi ) = ]^C 2 a + 0^ e ) (6) 

where <fi 2 is considered to be the smallest parameter in the problem. As a result, the cusp 
TBA system simplifies. Only the leading orders in the expansions of the other Y-functions 
are relevant. These are denoted by 
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The corresponding reduced TBA equations become [15] 
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and log F a = K a * log 
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(9) 
(10) 

(11) 
(12) 

(13) 



u=0 



7 The general near BPS case, with 6> 7^ has a similar degree of complexity. However, we leave that 
generalization to a future work. 
8 Namely, Y a<s (u) = F a ,_ s (-u). 

9 Here, F a is obtained from the one in [15] by changing the integration variable to u. 
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The hat on * is a convolution over the range —2g<u<2g. Our conventions for the kernels 
are 

*•<"> = ^CT) • *•<»> = {WSi KM ' RM = JW^ K ° (U) (14) 

and lf a & is given in appendix A. The cuts in K a and K a are chosen to be between —2g and 
2g. These integral equations are supplemented by boundary conditions. That is, at large 
real u the function Y m approaches its asymptotic value 

"-^('-^t) ■ (15) 

Finally, the energy is also dominated by the value of Y at0 at the double pole and reduces 

to 

1 00 C 

£ = -<P 2 B(\), B(A) = --£-==^_. (16) 

1 a=i V 1 + !% /« 

Note that even though the reduced TBA equations (8)-(12) are simpler than the general 
ones [15], still they are an infinite set of non-linear integral equations. Moreover, for L = 
the C a 's become oscillatory and truncating their sums becomes problematic for numerics. In 
the next section we will reduce that infinite set of equations to a finite set (FiNLIE) along 
the lines of [14]. As we will see, doing numerics with that finite set of equations is very 
simple, even for L = 0. 



3 Reduction to FiNLIE 

We will now reduce the Bremsstrahlung TBA equations (8) - (12) into a few functional 
equations for only three functions. We will mainly follow the methods developed in [14,23], 
based on Y-system or equivalently on the underlying integrable Hirota dynamics. We will 
start by applying the method to the infinite set of equations for Y m (10). 



3.1 Reduction of Infinite Set of Equations for y m 

We first consider the equation for the leading order of Y m , which are denoted as y s (13). 
The integral equation for them is 

log^m = Im,nS * log ( 1 + y J + <W-S* log — . (17) 

This is exactly the type of equation considered in [14,23,29]. That is, the integral equation 
(17) is equivalent to the functional equation 

log (y+y~) = Im,n log + 5 m , 2 log | (18) 
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plus the assumption that y m is analytic inside the strip |Imw| < 1/2. This is the Y-system 
for one dimensional wing (with the source log$/\l/ at its boundary). This equation can be 
solved in two steps. First, one introduces T-functions 

y m = J r+nT f l - 1 . (19) 

' m+l ' m—1 

The Y-system equation for the Y-functions is equivalent to the Hirota equation for the T- 
functions. The advantage of using the T-functions is that the general solution to the Hirota 
equation is known. That is, the functional equation (18) is automatically satisfied for m > 2 
provided that T s can be written in the form 



Ql(«+f) Q2(U + 

Q-2{u- l i) Qi(u- 



(20) 



for any four functions Qi,Q2,Q2,Ql- This change of variables thus already reduces the 
infinite tower of the 3^-functions to the set of four Q-functions. 

To reduce further the set of independent functions one can notice that the map (19) is 
not one-to-one. Indeed there is a "gauge" freedom in the choice of the % for a given y s . 
This freedom can be used to set two of the Q-functions to 1 

Q2 = Qi = 1 (21) 

Furthermore, due to reality of the Y-functions, we can always choose the remaining functions 
to be complex conjugate of each other 

Q x = -Q- 2 = Q and therefore T m = Q{u + f ) + Q(u - l f) . (22) 

This leaves us with essentially one single function! 

Next, we have to ensure the analytical properties that follows from the integral TBA 
equations (17). One of these is the analyticity of y m inside the strip — m f^ < Imu < 
This property implies that one can always choose the remaining function Q to be analytic 
everywhere above the real axis. As a result, we can write the following spectral representation 
for it 

00 

Q(u) = -iu- ! ^LJ^L for Imu > . (23) 

—00 

Here, the first term, —iu, follows from the asymptotic behavior of the Y-functions at large 
real u. The overall normalization is of course a gauge choice. 

Another analyticity constraint follows from the TBA equation for m — 2. It follows from 
that equation that y m should have at least two branch cuts starting at u — ±(i/2+2g) [13,14]. 
This is automatically true if p(v ) is nonzero inside the interval [— 2g, 2g] and vanishes outside. 
That is, the ends of the integration interval will create branch points in Q at u = ±2g. These 
cuts then translate through (20) into cuts at u = ±(i/2 + 2g) for 71 and hence for y 2 (19). 
In what follows, we will assume that the function p is real. As we will see, that assumption 
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is consistent with the equations (see also [14] for a proof related to Z 4 symmetry of the 
worldsheet action). 

To summarize, the infinite set of integral equations (17) can be packed into one single 
finite support function p(u) as 

T m = m + K m *p . (24) 

The function p is fixed from the functional equation (18) at m = 2 and the boundary 
conditions (15). In terms of p, that equation reads 



$ _ Ti(u + i/2 


f iO)Ti(u 


-i'/2-i0) 


(1 + $t 


*p- \p)i i +$i *p-\p) 


* Ti{u + i/2 


-i0)7i(u 


-i/2 + iO) 


(i + n 





(25) 



Here IfC 1 denotes the principal value integration with the kernel K s (u±i/2). The derivation 
of this equation is straight forward and is described in detail in [14]. 



3.2 Reduction of Infinite Set of Equations for X m 

Next, we move to the order <fr 2 piece of Y m , namely X m (13). At order <p 2 there are two new 
sources in the right hand side of (10) 

2 5 m>2 s* [$ - *] and 2 vrs(u)C m . (26) 

The first term is of the same form as the source in the equation for y m (17). It shifts the 
left hand side of (25) by 2 [$ — ty] and therefore is accounted for by a corresponding shift 
of the density p — > p + 4> 2 q. 

To relate the integral equation for y m (17) to the functional Y-system equation (18), we 
assumed that y m is analytic inside the strip |Imw| < 1/2. That is, for an analytic function 
f(u) inside the strip (including boundaries) we have 

f = S*g ^ f+ + f~ = g. (27) 
Suppose we modified the function f(u) — > f(u) so that it now has poles at u — ±i/2 

f(u) ee f{u) + - *W = f(u) + C«K x {v) . (28) 

Due to the contributions of these poles at the boundary of the strip, (27) is modified to 

f = s*g + -KCs{u) (29) 

where 

g{u) ee f + (u-i0) + f-(u + i0) = g + n CK 2 (u) + tt C 5{u) (30) 

The equation for Y m , (10) is exactly of that form! That is, the second term in (26) tells us 
that Y m should have poles at u = ±i/2 with residues ^i(p 2 C m /2 correspondingly. 
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Having understood this, our task becomes very simple - we have to modify our ansatz for 
Tm (24) to incorporate these poles. That can be easily archived by adding poles at u — |Z 
to the Q-function with to-be-fixed residues 



Q(u) = -i 
This procedure gives 

where 









i 













' dv p(v) + (j) 2 Q J_ b a 

2ni u — v 47r ^— ' it — i 



for 



lmw>0 (31) 



r 2 s 
\mu — m 



/12] + K m *g + 22 [bnK n - m (u) + b^ n K m+n (u)} (32) 



n=l 



Y, 



rrp "I (33) 

automatically solves the functional Y-system equation. Here, the 6 n 's are related to the 
residue of the poles. The polynomial — i[u + <p 2 u 3 /3] in Q is needed to support the boundary 
conditions at large u (15). 

It remains to fix the function g(u) and the residues b n . We fix the residues in the next 
section. The equation for g can be obtained by correcting the r.h.s. and l.h.s. in (25) and 
expanding to leading order in 2 . It reads 



rtp- {l + $t*P)Q , n p - (1 + $ l *p)g 



(34) 



Note that this equation also contains terms linear in T\. It is therefore a linear equation for 
g. 



3.2.1 Fixing Residues of Q 

We will now fix the coefficients b a . Following the construction above, logY m should be 
analytic in the strip |Imit| < 1/2 with only two poles at its boundary u = ±z/2. We have 10 



Y r . 



y m + 4> 2 y m Xm = y m + <p 2 {i + y„ 



T T 
in _|_ in 



' m+1 



T+ T - T 



m+1 



' m—l 
Tra—\ 



(35) 



The terms containing r+ and r m have potential poles at the origin. The condition that the 
residues cancel reads 

/i i — /i /i i — /i 

. (36) 



b m -i ~ ^m+i b m ^i — b m+ \ 



T m (i/2) r m (-i/2) 

That condition is automatically satisfied provided that the function T m is even (as the vacuum 
Y-functions are). 



5 All equations above should be understood to linear order in cf> 2 only. 
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Next, we proceed to the poles at u = ±z/2. The residues at these poles are fixed by (10) 
to be 

C 



logY m (w) 



2? 



u 



u + 



2 J 



(37) 



This has to be matched with the residue imposed by our ansatz. Close to u = i/2, Y m (35) 
diverges as 



Y m ( M )~ [i + y m (D] 



b m +2 — b m _^_b m — b m -2 



7m+l(§, 
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By comparing (38) with (37) we get a recurrence equation for the 6 m 's 
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This equation can be simplified a bit by introducing the numbers 

Co = p(0) , c m ^ = 7^(0) . 
As p(u) is also an even function, the real. In terms of these we have 



T m (ia/2) 



y m (ia/2) 



2 (c a — m+l C a _)_ m __i) (c a _ m _i C a _|_ m _|_i) 

In term of the c m 's, the recurrence relations (39) reads 
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and m = 2, 3, . . . . 



(38) 



(39) 



(40) 



(41) 



(42) 



This concludes our reduction of the infinite set of integral equations (10) for Y m {u) to a 
finite set of equations (25), (34), (42) for p(u), g{u) and for the coefficients b m . In the next 
section we will rewrite the remaining TBA equations by plugging our ansatz (32) into them. 



3.3 Simplified Equation for A m 

In this section we will analyze the equation for A m (12). Using the ansatz of the previ- 
ous section for Y m , we will get rid of the infinite sums in (12). That procedure is called 
"telescoping" in [14]. 

One of the terms in the equation for A m is K ab * log (1 + 34)- Using (19), those terms 
read 



K ab * io g (i + y b ) = jr K ab * [log r b + + log r b - log r h + x - io g 7^1] (43) 

b=2 b=2 

00 

J2[^ b -+K; b + -K a ^ 1 -K aM[ 

6=2 
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where we introduce a new notation for the shift f ±± = f(u±i/2±i0) so that in the first two 
terms the shifts are ±z/2 =F iO . The last two terms appear when we shift the summation 
index by 1. In fact K a ^\ = 0, so only one of these is non-zero. 

Now we notice something nice - the expression in the square brackets is almost zero. It 
is nonzero only due to the small displacements by ±i0 which may produce a delta function. 
Furthermore, one should be careful with the boundary case 6 = 2 where the branch cuts in 
K^2 approach the real axis and we get an extra contribution due to the square root cut 



Kl-iu) + K- + (u) - K a ,Au) - K M (u) 



8 a ,b8(u) + 6 bj2 x{u)K a (u) (44) 



where x{ u ) = 1 f° r — %9 < u < 2g and is zero otherwise. Thus the above integral simply 



gives 



jr K ab * log(l + y b ) = log T a - K a .2 * log Ti + KJ log T 2 . (45) 

6=2 

Let us rewrite the term K a p * log T\ by writing the kernel more explicitly and shifting the 
integration contour 

-K a ^\o Z T x = - l -(K+-K- + Kt + K-)*\ogT l (46) 

00+2O co — iO 

J du (R a + K^j log Tf + \ Jdu^Ka- K^j log 7?- 



2 

— oo+«0 — oo— iO 



Note that if the sign of the iO's was opposite, (46) was just zero. That is 

oo— iO oo+iO 

J du [k a + K^j log Tf = and J du [R a - K^j log T x + = (47) 



-oo— iO — oo+iO 



because the integration contour can be deformed to infinity as the integrant is analytical 
below/above the real axis. The difference between the two signs of iO is the contribution of 
the square root cut. That is, we have 

- K a , 2 * log 71 = \kJ log - l -K a i log {rr^-Ti-T^) (48) 

where again 

r +± (u) = T(u + i/2± iO) and T _± (u) = T(u - i/2 ± iO) . (49) 

In writing (48), we used the facts that the square root cut is taken between —2g and 2g. 
That is, 

K (u + iO) = { ~ ' 0) f ° r - 2 9<u< 2g 

ay ' \ +K a (u - iO) for it < — 2g || 2g < u 
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Now, the first term on the r.h.s. of (48) is nothing but ^K a *$>/ty according to (25). It 
cancels precisely with the term in (12)! That is (12) reduces to 



A, 



i~ *$7? r a 

-K a * log z t — h log — 

2 ^r 1 - + T 1 + -T 1 -r 1 ++ a 



(50) 



u=0 



We will see later that it is highly advantageous to introduce a new notation 

^T 2 §T2 



V 



'l 'l 'l 'i 



(51) 



where the second equality in nothing but (25). Furthermore, utilizing that 7^(0) = c a we 
have 



A a = K a *\ogr] + \og 



(52) 



That is a much nicer equation with no infinite sums! In the next section we summarize the 
FiNLIE and discuss its numerical solutions. 



3.4 Equations for Fermions and <£>) 

Our goal is to have a closed set of equations for a minimal set of functions. We will now 
exchange the two "fermionic" functions 11 \1/ and $ for the single function i] (51). As we 
discussed in the previous section there are two equivalent definitions of r\ (due to (25)) 



n 



{+> 



where we rewrote (51) in terms of the principal value 

2s 



{+> 



{-} 



(53) 



T} ±} = l - {r x ±+ + 7*-) = 1 + jdv p(v) K 1 (u-v± 



(54) 



-2fl 



The functions \1/ and $ enter into the equations through their sum and difference. We can 
now use (53) to express them in terms of p and r\ 



$ _ $ 



^ + $ 



\p 2 + 2T 1 { - } T 1 {+} 



{+} 



(55) 



What we have achieved by that is a set of closed equations for only three functions {p, r], g} 
and one infinite set of real numbers {C a }, as we will now summarize. 



ii 



The name is purely historical and does not refer to Grassman variables. 
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3.5 Summary of Bremsstrahlung FiNLIE 



In this section we summarize our finite set of nonlinear integral equations (FiNLIE) for the 
functions rj, p and g. As we will see, only the value of these functions on the interval [— 2g, 2g] 
will matter. This is particularly neat for the numerical implementation, which we discuss in 
the next section. 

The finite set of non-linear equations for rj, p and g reads 



1 00 



0=1 



ip 2 + 2T 1 { - } T 1 {+} 



x s * 



-2- 



+ n(K+-K-)C a -n6{u) d 



-{+} /V 



g = p- 



{+}< 



1 1 si 
4 1 4 



{-}2 



{+}' 



(7?" 



r. {+} ) ( p i 



T i-} T { + } 



2+2L 



-l) a ac a Jl + 



16g 2 Ag 



exp 



K a * log 77 



sinh(27ru) 



2nu 



(F2) 



:f3) 



:f4) 



Here, (Fl) and (F2) are the fermionic equations (8) and (9) written in terms of p and r\. 
The combination X 2 j (1 + 3^2) is given by (35) 



%T X 



1 + 3^2 V + Tf-\T 3 Ti 77 T 2 



r 3 , Tl TV 



2 T 2 



(56) 



Equation (F3) for £ is obtained from (34) using (55). Finally equation (F4) for the infinite 
set of coefficient C a is (11) where we used (52) for A. Note that the function 7] has combined 
neatly with F a from the boundary dressing phase (13). 

These equations allow us to find our three main functions {p, g, 77} in terms of T m arid 
r m . These functions themselves can be computed from p, g and C a through (24) and (32) 



T m = m + K m ip 



m 

T m = + mU 



(57) 



where the coefficients c m and b m are determined by (40) and (42) 



r m (o) 



b a +2 - b a 



2 2\ I n=1 

[ C a+2 ~ C a) X 



E 7rC2n-l _i_ 7rC2n 
„ 4c 2 „_i ^ 2-j 4c 2 „ ' 



n=a/2+l 



a/2-1/2 



- E 



7rC2r7 



n=0 



4C2n + l 



a G 2Z 
, a e 2Z + 1 



(5f 
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Note that we only need T\, t 2 and r 3 . In order to compute these, it is enough to know the 
difference b a+ 2 — b a : 

j oo 

n = -— +u 2 +Ki*q-) j (b n+2 - b n ) K n+ i , t 2 =t 1 + +t 1 ~ , r 3 = r^+^+r" . (59) 

n=0 

In the next section we discuss how this system of equations can be implemented numer- 
ically and present our numerical results. 

3.6 Numerical Solution of FiNLIE 

The numerical implementation of the Bremsstrahlung FiNLIE system, summarized in the 
previous section, is much simpler than the one for the full spectral FiNLIE [14, 23] . There 
are several reasons for that. First, there are only three functions (p, g, rj), that have to be 
updated by iterations. Second, only their values on the cut [—2g, 2g] have to be stored. For 
example, to efficiently store these functions at each step of iterations one can use a basis 
of Legendre polynomials, orthogonal on the cut. This makes it really easy to write a code 
which iterates these equations. 

There is however one new complication that appears in the system. The complication 
only comes about for L = and is due to the infinite set of coefficients C ffl . For L = these 
coefficients do not go to zero for large a. As a result, a naive truncation of the equations at 
some large a would lead to a wrong result. To overcome this difficulty we have computed 
the large a expansion of the C a 's in terms of rj analytically, up to order 1/a 20 . Having an 
analytic expression for their expansion at large a allows us to compute the tails of the sums 
over C a analytically. The remaining contribution can be truncated efficiently with essentially 
arbitrary precision. For the evaluation of the large a coefficients, we use (F4) which gives 
the C a 's in terms of some simple integrals of r\ and p. 




-1,0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 -1.0 -0.5 0.0 0.5 1.0 



Figure 2: Plots of rj,p and g at L = for various values of coupling from g = 0.05 (blue) to 
g = 0.25 (red) with a step Ag = 0.02. 

Since solving the TBA equations numerically is a rather standard procedure we describe 
our numerical algorithm very briefly. As a base for the iteration process we took the asymp- 
totic values of our functions po = £?o = and go = 1. At the n'th iteration step, we first used 
(F4) to reconstruct several coefficients in the large a expansion of C a . Then we plug the 
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Figure 3: Plots of rj,p and g for various values of coupling from g = 0.06 (blue) to g = 0.5 
(red) with a step 0.04 with L — 1. 

B(g),L=0, 1,2,3 

0.00 
-0.02 
-0.04 
-0.06 
-0.08 
-0.10 
-0.12 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Figure 4: Plots of energies for various L's as a function of coupling. The dots (from bottom 
to top) correspond to L = (blue), L = 1 (red), L = 2 (yellow), L = 3 (green). 

current functions p n , g n , rj n into the r.h.s. of (Fl), (F2) and (F3) to get the next iteration 
Pn+i, Qn+i, Vn+i- The accuracy of 10~ 7 is reached typically after 20 iterations and takes 
about 10 min on an average modern computer. 

The output of the iteration process are the functions r], p, g which are plotted in figure 2 
for L = and in figure 3 for L = 1. The goal of the iteration process is of course the 
energy. It is given in terms of the C a 's through (16). In figure 4 we presented the result 
for L = 0, 1, 2, 3 for a wide range of couplings. In that figure, we also plotted our analytical 
results (solid lines) which are derived in the next section. 

4 Analytic Thermodynamic Bethe Ansatz 

In this section we solve the FiNLIE analytically. We do so by understanding the analytical 
properties of our basic quantities rj and p. 
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4.1 Analytic Ansatz 



In the previous section we formulate a closed system of equations defined on the interval 
u G [—2g, 2g]. It turns out to be very useful to think more broadly and analyze the analytical 
properties of these functions on the whole complex plane, where the power of complex analysis 
can come to rescue. 

Consider first equation (Fl) 

oo 

pv = -7iJ2K a C a . (60) 

0=1 

It tells us that the combination p rj has only two branch points at ±2g and numerous poles 
outside the real axis at ia/2, a 6 Z. What we would like to know is what are the analytical 
properties of p and rj separately 

Unfortunately, it seems that there is no simple way to answer this question. At this point 
we can only argue at the emotional level that it would be very awkward if both 77 and p 
would have some other branch points that annihilate precisely when we multiply these two 
functions. Note that p must have branch points at ±2g as being a discontinuity of T (this is 
also clear from (F2) that rj does not have a cut on the real axis). Let us then be optimistic 
and assume the simplest possibility. Namely, we assume that 77 is a meromorphic function 
on the whole complex plane. We justify this assumption a posteriori by consistency with the 
other equations and by agreement with numerics. From (F2) it follows that 77 — > 1 at large 
u. Therefore, more explicitly, we assume that 

-i=n— < 6i > 

a 

for some set of roots {u a } and poles {v a }. In the next section we show that v a = ia/2. These 
poles of rj are responsible for the poles on the right hand side of (60). 



4.2 Fixing Poles and Residues of 77 



We fix the poles and residues of 77. We start from equation (F2) written in the form 



v- 



r i-} T i+} 



4 



r 2 /2 



s * 



i+y 2 



r/ + p + + rj p — 5(u)irC-[ 



Using (27), we convert it into the functional equation 



f)(T+ + ^^^ 



2 ■ 



+ c.c. = -2- 



X* 



By plugging in the explicit expression for (56), on the l.h.s and using that 
7i + = on the r.h.s, we arrive at 12 



T 2 + T 2 - 

12 That relation follows from (22). 



V 



+r — 



T3 n 



V 



(62) 



(63) 



(64) 
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This form is very useful for analyzing the poles of rj. First, note that T m — m and 
therefore T m can not have zeros at sufficiently small g. u As the function rj is analytic in g, 
the only potential poles in (64) are the following - poles of t s and {r]p) ± on the r.h.s and 
poles of Tj^~ on the l.h.s. Both type of terms on the r.h.s can only have poles at ia/2 for some 
integer a. It follows that also r] has poles at ia/2. However, this does not mean that all the 
poles of rj are at ia/2 since there could be another infinite series of poles which cancel in the 
finite difference. This possibility is excluded below. At the moment let us write a recursive 
relation for the residues of rj at ia/2. We denote them by e a 



Res rj = . (65) 

u=ia/2 2lX% 



Using that 



7rC a ba—s ^o+s I ■ /r>\ C a _|_ m Ca—m 

Res pr] = -7— , Res r s = — , T m [%a 2) = (66) 

u=ia/2 2m u=ia/2 2m 2 



we get 



(c a -l — C a _ 3 ) (c a _ie a _i + 7T C a _i) (c a+ i — C a+ 3) (c a+ ie a+ i + 7T C a+ i) 

H r = . (67) 



C a -1 (c a +l — Q1-3) c a+l (c a -l — Q1+3) 

It is very appealing to solve this second order recurrence equation by setting e a = —n C a /c a . 
However, we have to exclude other possible solutions by fixing the initial conditions. 

Consider first (67) at a = 1. We have Co = as the sum in (60) starts from a = 1. 
Moreover, eo = 0. Otherwise, it would follow from (60) that p(0) = 0, which generically is 
not true. Hence, the only solution is 62 = ~ E ^- As a consequence, for any even a we have 

p — nC a 

°a ,. 

It is left to fix the residues {e a } at odd a's. The first of these, for a = ±1, are at the 
boundary of the strip —1 < Im-u < 1. Therefore, we go back to the integral equation (F2). 
It gives us direct accesses to the analytical properties of rj inside this strip. The only source 
for the poles in the r.h.s. of (F2) is the term with 5-function (which so far played no role). 
It gives 

rrr+ + g ^ _ ttq 

V T2/2 ~ 2cosh(vra) 1 ' 

comparing residues at i/2 of right and left hand sides we get 

r . r . c 3-ci _ ECl( r _. _ C3-C1 N si 

0161 2 2 {Cl = — t\C\ =► Cl = . (69) 

Cl 



C3+C1 

4 



We conclude that for any a 



(70) 



13 The function T m only depends on p in the interval [— 2g, 2g]. Moreover, at g — > the solution reduces 
to the asymptotic one (15) where rj ~ 1, p ~ and therefore T m ~ m. 
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Note that (F2) implies that there are no other poles in the strip (F2). It follows from (64) 
that there are no other poles in 77 except ia/2. Furthermore, taking into account that form 
(F2) i] 1 for large u, we can write the following representation of rj 



V(u) = l-J2—K a (u) 



(71) 



This is the key equation which will allow us to solve the system analytically. This equation 
was successfully tested numerically! In the next sections we also constrain the zeros of f](u), 
thus writing a closed equation for rj(u) alone. 



4.3 Effective Baxter Equation for Zeros of 77 

In the previous section we have found that the poles of r](u) are at ia/2. Then (61) becomes 



■q(u) = ] I 



OO 9 9 

u - ul 



\u 2 + k 2 /A 



(72) 



where we also use that rj is an even function, so its zeros and poles should generically come 
in pairs. 

The goal of this section is to find the zeroes Uk in closed terms. For that we plug the 
ansatz (72) into (F4). We see that we can evaluate the integral in (F4) explicitly in terms 
of Uk- Using 



sinh(27ru) 
2nu 



n 

k=l 



u 2 + fc 2 /4 
k 2 /A 



(73) 



the exponent of (F4) becomes 



00 1 



k=l 



29 

2 2 

, ~ , u-u% 

-29 



(74) 



where we have written the convolution as a contour integral around the cut [—2g, 2g]. This 
allows us to blow the contour to infinity. There are two type of obstacles on that way - poles 
of K a at u = ±m/2 and logarithmic cuts starting at u = ±Uk- By picking the residues at 
the poles and integrating the discontinuities along the cuts, we obtain 



poles 



cuts 



log 



-a 2 /A - ul 



k 2 /4 



r](u) smh(27ra) 
hm log = log 

u-nii/2 2lXU 



c 

a 

ac a 



k=l 



.1 j 



(75) 
(76) 
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where the contribution of the cuts is written in terms of the Zhukovsky variable Xk = x(uk) 
and y a = x(ia/2). 14 Combining both contributions together we get from (F4) 



2+2L 



-l) a ac a L/1 + 



16c/ 2 Ag 



oo J2 1 



:-ir— U-^A 



(77) 



The expression in the large round brackets is just i/y a and (77) simplifies to 



2+2L 



oo „2 1 



n 



k=l X k Va 



(7E 



This strongly reminds some kind of Bethe ansatz equation! To solve it we use the technology 
of Baxter equation. Namely, we construct a function 



T(x) = a; L+1 Q0r) + ±JLq(1/ x ) , Q(x) = J] 



OO o o 

Uj J,, Uj 



k=l 



X i 



This function has very simple properties. These are 

T(x) = (-l) L T(l/x) , T(y a ) = 



(79) 



*0) 



and the absence of singularities, apart from essential singularities at x = and x = oo. To 
analyze the behavior at infinity (and at zero) let us relate the Q-functions with rj: 



Q(x)Q(l/x) = ]( 



00 2 2 

Wt - U 



,=1 ^ 



. sinh(27Tii) ~ 



where the constant C = YYk=i zj^ ■ When y — > ±00 we also have u — > ±00, thus rj(u) 



9 x 



Q(l/y) ~ 1, see (72) and (79). Therefore, we have 



ZTYU 



This allows us to estimate the behavior of T 



T(y) ~ y L sinh(2iiu) 



C 

2^g 



at 



y — > ±00 . 



The function which enjoys this asymptotics and satisfies (80) is 



C 



l+L 



T(y) = smh(27Tu)P(y) with P(y) = C lV L + C 2 y L - 2 + ■■■+ , 

yL 

where Ci+2-i = (— l) L Ci- Note that C = 2ixgC\ and therefore 



(82) 



(83) 



54) 



u Q(z)Q(l/s) 
(7 Ci sinh(27TM a; ) 



»5) 



4 The Zhukovsky map is defined as ^ = x + \ and we consider the solution outside the unit circle. 

g x 
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Finally, recall that that Q(x) = [Pl(x) sinh(27ru)/x L+1 ] + , where the subscript (+) stands 
for the regular part of the Laurent expansion. We can therefore write an integral represen- 
tation for Q in terms of Pl as 



/ dx P L (x) smh(27iu x ) 1 

Q(y) = - <p 2 — 2 1 ■ 86 

' 2m x l — y l x L 



\A>\v\ 



To fix T(x) completely, it remains to find the coefficients Cj. That is done in the next 
section. 



4.4 Solution of Effective Baxter Equation 

For simplicity, we will first demonstrate the procedure at L = 0. Then, we will move to the 
general case. 



4.4.1 Q-function for L = 

For L = the general equation (84) becomes 

T(y) = d sinh(27rw 2/ ) . (87) 

It is useful to expand T in Laurent series around y = 0. For that we use the well known fact 
that sinh is a generating function for Bessel functions 

oo 

T(y) = d (4*0) y 2r+1 (88) 

r=— oo 

comparing with (79) we see that the first term contains only negatives powers of y whereas 
the second term has only positive powers. That is, we have 

QM - ft ^ = CiE hr +1 (4«g) y* = - Cl I . (89) 

k=l k r=0 \ \t\ \ 

The constant C\ is fixed from the condition Q(0) = 1. It reads 

Ci = hh) ■ m 

We will now generalize the construction of Q to any L. 
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4.4.2 Q-function for General L 

For general L we had 



T(y) = y L+1 Q(y) 



^Q(l/y) = smh(27T Uy )P(y) 



(91) 



and therefore T{x) also has the zeros of P{y) in addition to y a = x(ia/2). The coefficients 
in the Laurent expansion of P(y), C r (84) should be determined so that T(y) has a gap in 
its Laurent expansion from — L to L 



T(y) = sinh(2vra) ]T C r+1 y L ~ 2r = ■■■+ + + ■ ■ ■ + x y L ~ l + y L+ 



r=0 



J 

/r -r< yl 



(92) 



More explicitly, the coefficients Cj satisfy 

L+l 



dx P(x) sinh(27ra a 



2ni x 



L+2-2i 



i=i 



'-l)% L+ i for i = !,...,£+ 1. (93) 



These are just L+l linear equations for L+l variables, the C a 's. We have C a = (— 1) L [A / 1 L 1 ] a L+1 
where the (L + 1) x [L + 1) matrix M. L is given by 



M< 



( I-i h 

1-3 1-1 

h-2L h-2L 

\I-l-2L h-2L 



The function P^ can be now written as 



I2L-3 

I-l h 

1-3 I-l J 



(94) 




(95) 



For example, for L = 0,1,2 we get 

x — 1/x 



Pa 



1 



Pi 



Pi 



hx 2 -(h + h) + h/x 2 

il - hh + hh - H 



(96) 



In terms of z = | [x + 1/a;), the function Pl(x(z)) is a polynomial of degree L for even L. 
Similarly, for odd L the function Pl(x(z))/[x(z) — l/x(z)} is a polynomial of degree L — 1. 
In appendix C, we show that these are orthogonal polynomials for some simple measure. 
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Having explicitly the function P^ and therefore also Q through (86), we can now easily 
reconstruct rj(u) itself via (85). In the next section we use the solution for rj(u) to compute 
the energy. In that computation, we will also need a function with one less degree than P^. 
It is defined as qiiy) = [P(y) sinh(27rw)] + — P(y) [sinh(27ru)] + . It follows that qiil/x) = 

Ui +l Q{yi)i see appendix B for details. 



(— 1) L 1 9i,(x) and at the roots of Pl we have q{y^ 
We find that 



/x 1 ~ L \ 



Ql(x) 



„3-L 



„L-1 
L+l 



\x L+1 J 








h 








hL-5 






l2L-l\ 



h 





/ 



(Cl+i\ 
C l 

V Cl J 



(97) 



Notice that what enters in the expression above is the part of Aii above the diagonal. 



4.5 Expressions for Energy 

We would like to compute the energy (16), repeated here for convenience 



Bjg) = --y . C (9f 



For that, we need to find all C 's. That is done through the function p(u). Indeed, note that 
by rewriting (60) using (71) we have 



= ^ a *ya~aW = _ (gg) 



and therefore 

B(g) = - lim u p(u) . (100) 

Moreover, from (99) we see that 

p(ia/2) = c a . (101) 

Combining this with the definition of c a (40) (and using that p is an even function), we arrive 
at 

P^IA = a _ + J*L pVL , aeZ +. (102) 
2 2 J 2iriu-ia/2 ' v ; 

-2s 

To better understand what this equation tells us we introduce the function 

H(u) = u+ Jdvj^l.m. m 
K J J 2-kv-u 2 y ' 

This function has very familiar analytical properties. First, due to (102) it has zeroes at 
v = ia/2. Second, it also has poles at the zeroes of 77, {uk}, due to the ip(u)/2 piece. Finally, 
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by its definition it seems like H(u) has a branch cut at the interval [— 2g, 2g] due to the 
singularity of the integrand. However, the discontinuity of the integral, ip, precisely cancels 
with the discontinuity in —ip/2\ Thus H(v) is a meromorphic function with poles at zeros 
of rj and zeros at ia/2 (which are poles of 77). As it behaves as ~ u at infinity we conclude 
that 

11 

H(u) = -F-r . (104) 



T}(u) 

We can now reconstruct p(u) from i](u). For —2g < u < 2g we can rewrite (103) as 

2.9 



U 



r){u) 



u 



dv p(v) 
2tt v — u 



(105) 



-2s 



Finding p is now reduced to a standard Riemann-Hilbert problem. Using that H(v) — v 
decays at infinity, we get for u inside the cut 



2S 



p{u) = 4 



-2s 



dv y 4g 2 — w 2 1 

27T ^4^,2 _ y 2 



U — V 



_T)(v] 



— V 



-2g < u < 2g 



(106) 



Analytical continuation outside the cut is 

2s 

p(u) = 



2dv\Ju 2 — Ag 2 v 1 — 77(1;) . 1—77(1/.) 



7T7 



-2a 



U — V 7](v) 



(107) 



Using (100), we find for the energy 



2a 



Bl(9) 



dv 



7T 



-2a 



V 7 ^ 



1-77(7;) 2 1 - ? ?H 

— hm m 



Tj(v) 



Tj(u) 



(108) 



At this point we can state the complete analytical solution of the initial problem of finding 
Bi(g) for any given L. In the next section we use this general solution to find a simple 
explicit expression for B^g). 



4.5.1 Closed Expression for Energy 

Let us begin by evaluating the second term in (108). For that we have to expand i](u) at 
large u, where x u is also large. Using (91) and (85) we have 

, \ u [sinh(27ra)P L (x) - {-lfx-^G&l/x)] Q(l/x) uP L (x) Q(l/x) 
Tj{U) = — 



g C\ x L+1 sinh(27T7i) gCiX L+1 
where the second step we drop the exponentially small terms. Next using that at large x 
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where 

L+l 

/j^Cj/^ (109) 

are the coefficients in the Laurent expansion of P(x) sinh(27ra). Putting the pieces together, 
we find 

,(«) = i + 1 (i + /„ +!)=>- ^y-^ = + /« + 1) • ("o) 

We now turn to the integral in (108). We write it as a contour integral 

(if f 2 1 — ri(v) f dx ( 1 \ u 9 , 

Tw-v = - 9 /^( 1 + ^)^)- 2g • (m) 

-23 

In this form it is rather pointless to try to compute this integral by poles since there are 
infinitely many poles at u = Uk situated roughly along the imaginary axis. Instead, the idea 
is to use again (91) and (85) to factorize this function into two parts analytical inside/outside 
the unit circle 

^ + (-1)V7. . (112) 



rj(u) x L+1 Q(x)P L (x) l/x L+1 Q(l/x)P L (x) 

Note first that each of these terms separately has poles due to the zeros of Pl{x) in the 
denominator. These poles, {yi}, are localized on the unit circle (see appendix C for a proof). 
They cancel between the two terms, so we can choose the integration contour in (111) to 
go slightly outside the unit circle. Second, at least at small enough coupling, Q(x) does not 
have zeros inside the unit circle. That can be seen by noting that at weak coupling rj — > 1 
and therefore, as follows from (85), the zeros of Q(x) can only be near u = ia/2. Other 
values of the coupling are defined by analytic continuation. Now, the second term in (112), 
analytical outside the unit circle contributes only by a residue at infinity. For the first term 
we contract the contour to the origin, picking the residue of the pole at the origin and the 
residues at the zeros of Pi. All in all, the integral in the r.h.s. of (111) gives 



-9 



dx f 1 \ u f . Co \ 9 ( 1 A Ci , 



2ni \ x 2 J rj(u) \ C x ) ^ \ yf J q(y i )d y P L (y l 



where we used that q{yi) = y^ +1 Q(yi), see discussion above (97) and appendix B. Combining 
all terms together, we arrive at 



(114) 



In the next section we will further simplify this expression. We will write it in a form 
which does not require finding the roots of Pl- Before doing so, let us consider two simple 
cases where L = and L = 1. 
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• For L = this equations simplifies considerably since Po(x) = is just a constant, 
so the sum over zeros vanishes. Moreover, C2 = and thus 

B L=0 = 9 \1 - fs) = 9 2 (l ~ hC t ) = g 2 (l - (115) 

in agreement with the localization prediction (3)! 

• For L > there is no result available from the side of the localization so it is interesting 
to see what our equation gives. For L = 1 things are still extremely simple. P\{x) = 
Ttg^j^- furthermore using that 

° 2 - 1 f - 9 /4 

we find 

Bl=i 

We see that the result for L = and L - 
Bessel functions. This holds true for any L, which is completely nontrivial from the expression 
(114). This indicates that (114) can be further simplified. 



2^ -2 „(„.\B.p T („.\ =2 —r- ( 116 ) 



yf q(yi)d y PL(yi) h 

1 n - hh 



(117) 



1 takes the form of a rational function of modified 



4.5.2 Explicit Expression for Energy 

We will now simplify our result for the energy (114) into an explicit combination of Bessel 
functions. 

First, note that our general expression for has the form 

B L (g) = g 2 A L (g)+g 2 r L (g) (118) 

with 



Al = 1 - /l,l+3 - tt^ and r i = y^( 1 + -2") 
Cl,i Vv ViJ 



qL(yi)d y P L (yi) 



(119) 



in this section we add an extra index L to Cj and to indicate that they are different for 
different L's. 

The key observation is that the complicated part tl can be expressed in terms of the 
simple part A, 

L-l L-l 

r L = 2j2(-V L ~ l A t B L = g 2 A L + 2g 2 Y,{-^) L ~ i A l . (120) 

i=0 i=0 

We observed this relation from the explicit expressions for Al and tl for first several L's. 
We did not find any simple analytical proof of this relation for general L, however we verified 
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the relation numerically with 200 digits precision for L = 1, . . . , 100 for several values of g. 
This relation allows us to overcome the step involving finding zeros of Pl{v) which is rather 
hard to do analytically for L > 3. Furthermore it allows to further simplify our result. Next 
we use that 



9i 

Ci 



detM^ 



detM 



L,L+3 



detA4+ 2 i 



detM 



(2,1) 
L+2 



detM 



(i.i) 

L+l 



detM 



(i-i) 

L+2 



121) 



where M^ is the matrix obtained by deleting the a th row and 6 th column of Ml- The first 
equality is clear from (95). The second equality in the relation for /l,l+3 holds true for any 
matrix whose entries depend only on 1 2i — 2j + 1 1 . As a result we get 



A; 



detM 



(2,1) 
i+2 



detM 



(i,i) 

i+2 



+ 



detM 



(2,1) 



detM 



(i,i) 



(122) 



By plugging the A^s into (120) we obtain our final result for the energy 



B L (g) = g 2 



detM 



(2,1) 



detM 



L-2 +2 det -^L 



(2,1) 
+1 



detM 



(2,1) ~ 



(i,i) 

L+2 



detA<£u 



detM 



(i,i) 



(123) 



which can be also written in the form (4). 



5 Classical Limit 

In this section we make an important test of our construction. We consider the classical 
limit L ~ g — > oo of our result (123) and compare it to the classical string prediction. In 
the large L limit our explicit result (123) becomes rather complicated since the size of the 
matrix Ml goes to infinity. We found it useful to rewrite (123) in terms of a matrix model 
type of integral to perform this task. We also discuss the algebraic curve, arising both in 
the classical open string and in the matrix model descriptions. We start from the string 
description and then move to the classical limit of our result. 



5.1 Classical Strings 

The classical string solution for general L,0 and 9 is considered in details in appendix E, 
where we derive its classical energy in a parametric form. In particular for small angles it 
gives 

L = Ag [K(oo 2 ) - E(oo 2 )] (124) 

E = L+ ' j(92 -^^m (125 > 
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and therefore 15 

The energy (126) is only one charge in an infinite set of conserved charges. All these 
charges are nicely encoded in the so-called algebraic curve. The algebraic curve is represented 
by the quasi-momentum p(x) - a function of the spectral parameter x that encode all the 
conserve charges in its Taylor expansion [30]. In general, finding that function at strong 
coupling involves solving a linear system of differential equations. Here we make a shortcut 
by relating our solution to su(2) Folded String, for which the algebraic curve is known [31,32]. 

5.1.1 Relation to su(2) Folded String 

The S 3 part of the open string solution dual to the Wilson loop is analogous to the su(2) 
folded string. To see the relation let us compare the global charges. For the folded string, 
these are the energy and two angular momenta on S 3 



£ folded = 4^i? 2 + -^-2cos(20) K(sin 2 (0)) (127) 
J™ = 4s (r + I) [K (sin 2 (0)) - E (sin 2 (0))] (128) 

4° Wcd = 4^i?-^E(sin 2 (0)) . (129) 

The result for the energy is given in a parametric form in terms of two parameters, R and 
defined by the equations for (J[ olded , J2° lded ). These parameters have a very clear meaning in 
terms of the algebraic curve. The algebraic curve has two symmetric cuts with the branch 
points at ±i?e v. We immediately see that the expressions for Jf° lded (128) and L (124) are 
very similar under the identification 

oj = sin0 . (130) 

In order to compare the folded string energy with the cusp solution we have to extract 
the S 3 contribution to the energy of the cusp solution in Appendix E. This is given by 16 

E s s =4gco K(co 2 ) + 0(6 2 ) . (131) 

We see that again the pattern of the elliptic functions is exactly the same and we have to relate 



the {^goS\ = [4p sin <p\ pre-factor to AgyR 2 + — 2cos(20) that looks more complicated. 

However, if we set R — 1 the square root simplifies to [8g sin 0] and therefore, up to a factor 
of 2, is exactly what we need. Similar mismatch by 2 appears in 

jfolded ( and jfolded = q) 

This factor of 2 has a clear origin - it is due to the fact that the folded string is a closed string 



15 In appendix F the small L expansion of this result is compared with fixed L, large g expansion of the 
TBA prediction. Despite the potential order of limits issues the results agree perfectly. 



16 Which is 2<7 s in the notations of appendix E. 
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solution made of two folds, whereas the solution dual to the Wilson loop is an open string 
solution associated to a single fold. Having matched the parameters of the two solutions, we 
immediately read the S 3 part of the algebraic curve 17 . 

5.1.2 Algebraic Curve 




Figure 5: Algebraic curve ^—^p{x) (real part is plotted). The cuts are discretized by zeros of 
Pl(x) (red dots). 

Above we identified uniquely the (S 3 part of) algebraic curve corresponding to our clas- 
sical solution. It should have two cuts on the unit circle with the branch points at ±e ±w ^ 



where 



arcsinur . This algebraic curve was studied in details in the relations to Giant 



Magnons, which are a particular limit of the folded string solution. The quasimomentum is 
known explicitly and for \x\ > 1 it reads 



p(x) 



7T + 2K 



-x 2 + e 2 ^ 



sin 



sin 




x 2 e 2i<t> / 2ix sin 



ar 



1 



1 ) + -^-F 1 -4cos0KF 2 (132) 
cos 6 



e 2i<f> _|_ 1) (g 



X 



2i<f> 



1) 



+ X) 



2i<b 



icf> 



X 



(e 2 ^ - 1) (e^ + x) ' 



tan 



tan 



where K = K(sin 2 0), E = E( sm 2 (f>). For \x\ < 1 one should use the relation p(x) = p(l/x). 

One can see that the function p(x) indeed has two cuts (see figure 5). Usually the 
quantization of the classical limit consists of a discretization of the cuts by poles. For 



17 The AdSs part should be identical to the S 3 part in order for their contributions to cancel in the BPS 
limit. 

18 Curiously the S 3 part is precisely the Giant Magnon with momenta 4>. 
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example, in the spectral problem the cuts are collections of the Bethe roots which become 
dense in the classical limit. In our case of the cusp anomalous dimension, it seems naively 
that there are no Bethe roots at all since we are dealing with Vacuum - the reference state 
on top of which the traditional Bethe roots correspond to excitations. From this point of 
view the situation is very different. However, in the next section we argue that the cuts 
are formed by roots of the polynomials Pl. These are the holes yi in the effective Baxter 
equation (79). 

Our interpretation of this is that the standard asymptotic Bethe ansatz at the cusp is in 
a sense an effective Bethe ansatz build over a nontrivial Vacuum. That vacuum, by itself, is 
quantized and is characterized by an additional set of roots t/i (zeros of Pl). These "hidden" 
variables yi are not visible in the standard ABA equations. Their existence is represented 
by the boundary scattering matrix and their contribution to the energy only comes about at 
finite size corrections. As a result, in the standard ABA equations that information is wiped 
out. In particular the ABA energy is exactly zero. 

Another important difference between the spectral problem and the situation here is the 
following. Here, even in the classical scaling limit, the standard ABA equations do not de- 
scribe the system exactly. This is contrary to the spectral problem where the Bethe equations 
describe accurately the classical limit including all finite size effects. That difference should 
disappear when the extra level of yi is included into the Bethe equations. In particular, as we 
will see, the ground state energy is described accurately in this case. It would be interesting 
to write the full system of Bethe equations which include the y^s. In other words, it would be 
interesting to include the excitations on top of the vacuum into the effective Baxter equation. 

We will now turn to the classical limit of our result (123). 



In this section we give an interpretation of the algebraic curve, discussed in the previous 
section, in terms of the quantities arising in the exact quantum TBA description. As we will 
see, the key object in the analytical solution of the TBA equations is the polynomial Pl{x). 
In the classical limit, the roots of Pl will condense and form the "vacuum" curve. 

5.2.1 Matrix Model Reformulation 

In order to compute the classical limit of our final result (123) it is very convenient to 
rewrite the determinants in a form of multiple integrals. For that we use the following 
integral representation of the modified Bessel function 



5.2 Classical Limit of TBA 




(133) 
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then, for example, det M L _i becomes 



det M l- 



n 



dxk 
2m 



-2irg{x k +l/x k ) 



.X j 
.X 9 



„2-2L 

'L-\ 

,-2L 



.X -y 
,X 9 



r 4-2L 

l L-l 

„2-2L 



2L-6 2L-4 

2L-8 2L-6 
x 2 x 2 



.X 2 



(134) 



x 



X 



L-l 

-2 



the determinant under the integral in the r.h.s. is related to a Vandermonde determinant in 
a simple way 



L L 

dXh . , ox -2KgJ2( x k+l/xk) 



det = i J] -_^A(z 2 ) e * (135) 

where 



A(.rj) ( 136 ) 

To bring the expression to the form of the partition function of a matrix model we have to 
get rid of the x 2 ^~ 2 in the denominator. For that we use that A is antisymmetric, so we can 
(anti)-symmetrize this factor to bring it to a Vandermonde determinant 

£(-i) M n^ M = !r!l (137) 

so that 19 

, /A dx k A 2 (X 2 ) ~27rgf:(x k +l/x k ) 

k K 



As it is common in the matrix models let us introduce the resolvent 

Gi(i) "d^Mi^/nss|r— i?^i e • (139) 

Note that det = detM L _i, so in order to get the generalized Bremsstrahlung function 

we need det -M^ 2 = det Ai^ . This can be also expressed in the form of the integral (135) 
with an extra \ factor which in terms of the resolvent becomes 



M 



(2,1) 
L 



M 



-LG' L (0) . (140) 



So far we used no approximation. In the next section we evaluate this quantity in the classical 
limit g ~ L — > oo. 

19 This ensemble also describes the 0(— 2) model in constant magnetic field near the critical point [33]. 
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5.2.2 Saddle Point and Algebraic Curve 

In the large L limit we have to find the saddle point for the integral (139). The saddle point 
equation is 

xl - 1 / 1 1 \ L . . 

At the saddle point where (150) is satisfied, the resolvent is simply given by 

Ci(a0=!G -( a0s ^g(_L- + _l_). ,142) 

Let us now show that in the classical L ~ <? — >■ oo limit, the solution to the saddle point 
equation is given by an algebraic curve with two cuts. In fact we will see shortly that this is 
precisely the classical algebraic curve arising in the classical string description discussed in 
the previous section. We introduce the function 

P(x) = -^-r - -^-Gf(x) (143) 
g x z — 1 g x z — 1 

in terms of this function the equation (150) simply becomes 

p(x k -HO) +p(x k - iO) = -2tt . (144) 

This is an equation familiar from the general classical algebraic curve approach [30]. Fur- 
thermore one can see from (150) that the set of roots {x{\ is invariant under inversion. This 
implies that 

p(x) = p{l/x) . (145) 
Together with the other asymptotic properties 

p(-x) = -p(x) , p(x) ~ -— + 0(l/x 3 ) (146) 

gx 

the function p(x) is uniquely determined. We can reconstruct the functions solely from these 
properties along the lines of [30], however all these properties are exactly the same as the 
properties of p(x) in (132)! So we conclude that p(x) we defined in terms of the roots in 
(143) is the one written in (132). This implies that we can immediately extract the resolvent 
and thus the ratio M^/M^ in (140) 



M 



(2,1) 
L 



g (x 2 — 1) p(x) L 
2x 2 2a; 



L 
3" 



2 E(w 2 )+K(^ 
1 + uo l 



E(u 2 ) - K(u 2 ) 



(147) 



where oj is an implicit function of ^ via K(w 2 ) — E(w 2 ) = ^ 

Now having this expression computed we obtain the generalized Bremsstrahlung function 
from (123) 

M (i,2) 

B d = -9 2 dl%v, (148) 



,2^2 

' M 



(1,1) 

L 
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using that d L = 4g ^ 2) ^ we get 



in the perfect agreement with (126)! 



5.2.3 Zeros of Pl as Quantization of Algebraic Curve 

In this section we will show that in the classical limit zeros of Pl(x) coincide with the solution 
to the saddle point equation (150). Therefore, they can be thought of as the quantization of 
the algebraic curve. 

Note first that up to a constant factor, Pl-i{x) is given by a determinant of Bessel 
functions and the monomials of x (95). Therefore, using (133) it takes the same form as the 
matrix-model type of the integral (134). The only difference is that in the last line in (134) 
we write x instead of xl and remove the integration over xl- In the classical limit the integral 
is saturated by the saddle point values Xi — > xf (150), so that Pl-i(x) becomes proportional 
to the determinant with fixed Xi = xf. Now it is obvious to see that the determinant is zero 
whenever the argument x is equal to ±xf . We conclude that in the roots of Pl(x), denoted 
by yi, satisfy the saddle point equation 

L 2 ^ 1 
g^L^oo : - V = n k ir g , k = 1,. . . ,2L (150) 

Vk - 1 Vk ~ 1 J^ k Vk - Vj 

where = sign Key k . It means that the zeros of Pl{x) discretize the curve and allows us 
to interpret them as the quantization of the classical curve. At the same time, these zeros 
are the "holes" in the effective Baxter equation (91). 

Note that (150) formally coincides with the strong coupling limit of su(2) spectral Bethe 
ansatz in the case where the number of roots is equal to the length (which is equal to 2L 
due to the open/closed double folding)! This seems to suggest that the initial vacuum with 
Z L at the cusp as being a maximally excited state of some other reference state. It would 
be interesting to understand this point in more details. 20 

Having established the classical limit of the zeros of Pi (150), we can now write an explicit 
expression for the classical limit of Pl itself. That is done in appendix D. 

An interesting question is whether there is a simple generalization of (150) for any L 
and g. For example, for g = and any L we find that Pl{x) is related to the Chebyshev 
polynomials (see appendix C). Their zeros are known explicitly y k = e t7Tk ^ L and they satisfy 
an equation similar to (150) 

g = 0, anyL : - V — ^- = , k = l,...,2L. (151) 



20 Note that the choice of treating Z L as the cusp vacuum was very convenient as it preserve supersymmetry 
in the BPS configuration. Other possible choice is to use a complex scalar made of the same scalars that the 
Wilson lines couple to. Away from the BPS configuration, supersymmetry is anyhow broken and that other 
choice preserve more (bosonic) symmetry. 
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Note that (151) is similar to the SU{2) case [34]. 



6 Summary and Discussion 

In this paper, we have computed analytically the quark-anti-quark potential on S 3 in the 
near BPS configuration and with L units of R-charge. Let us summarize the main steps in 
the derivation and identify potentially important future directions. 

We started from the vacuum Thermodynamic Bethe Ansatz (TBA) equations in the so- 
called "simplified" form listed in section 2. By using the methods developed for the spectral 
problem (known as FiNLIE) [14, 23] we reduced this infinite set of integral equations to just 
three equations for three functions p(u),r)(u) and g{u). Quite remarkably, only their values 
on the cut u G [— 2g, 2g] enter the equations, making them very convenient for numerics. By 
solving them numerically for various values of the coupling, we confirmed with high accuracy 
the localization result for L = and made numerical predictions for other values of L. Next, 
by using assumption about the simplicity of the analytical properties of rj, p and by analyzing 
the structure of their poles and zeroes, we derived an explicit equations for them (85,106). 
These equations agree nicely with the results of our numerics, verifying the assumption! 

Perhaps, even more excitingly we found that the zeros of i] (or equivalently the poles of 
p) satisfy some kind of Baxter equation (79). This Baxter equation has a rather curious form 
- instead of the familiar quantum shifts in the spectral parameter u by ±i/2, we observe a 
"crossing" type of shifts x — > 1/x. Curiously, the same shift emerged before in the strong 
coupling Thermodynamic Bubble Ansatz (TBuA) equations for classical minimal surfaces 
in AdS [35-37]. Our construction on the other hand is completely quantum. That is to say, 
it is valid at any value of the coupling. It would be very interesting to see if our approach 
provides a natural quantization of the minimal surfaces linking these two approaches. 

Another important object, emerging in our consideration, is the polynomial Pl{u). Zeroes 
of Pl{u) are the "holes", i.e. nontrivial zeroes of the Baxter function T(w). We have found 
that the polynomial Pl{u) is related to a matrix model resolvent. This interpretation was 
especially useful in the classical limit L ~ g — > oo where the matrix model integral is 
saturated by a saddle point. Very naturally, the algebraic curve describing the saddle point 
of the matrix model coincides precisely with the classical algebraic curve of the integrable 
worldsheet theory. Moreover, we have shown that the roots of Pl{u) coincide with the matrix 
model eigenvalues in the classical limit. This implies that these roots provide a quantization 
of the classical algebraic curve. 

There are several interesting future directions of study: 

• It would be interesting to understand the physical interpretation of the matrix model 
we obtained. Note that this matrix model is different from the ones usually arising from 
localization techniques where one is dealing with N c x N c matrices. Here instead, we are 
working in the planar limit from the very beginning and the matrix model we obtain 
is for L x L matrices, where L is the number of scalar Z fields inserted at the cusp. 
One possible interpretation is that the matrix model is related to a quantization of the 
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transfer matrix in the physical channel. Such picture connects neatly with the TBuA 
for minimal surfaces discussed above where the elements entering into the equations 
are the transfer matrices between edges of a polygon Wilson loop. 

• Whereas for L = there is already known result from localization (which we suc- 
cessfully reproduced), for L > our results are new. The localization techniques 
of [20,21,26] can be generalized to the case considered here. 21 Our results therefore 
provide for the first time a prediction from integrability for the localization techniques. 
Having an alternative derivation of the results may shed some light on the relations 
between the two very different methods. 

• It may be interesting to find a similar regime of the spectral TBA to the one considered 
here. For example, the slope function [38] in ABJM theory [39,40] is one of possible 
candidates. 

• In this paper we considered the near BPS configuration in the case where the angles, 
4> and 9 are small. It should not be too hard to generalize our result to the case where 
((f) 2 — 9 2 ) is small but <fi is finite. In particular, it would be nice to understand the 
generalization of our classical curve. 22 

• It would be interesting to compare our results at strong coupling with the string com- 
putation at one loop. From the TBA point of view this would require solving the 
matrix model to the next order (a part of the result is in appendix F). On the classical 
string side of the calculation, it should be possible to obtain the one loop result from 
the algebraic curve [41]. As we already know the algebraic curve, the result to a large 
extent can be read from [32,42,43] (up to a subtlety of open vs. closed strings). 

• We hope that the methods developed here could be extended to any value of the angles 
0, 9 and in particular to the general spectral TBA. The analog of the function 77 with 
simple analytic properties should also exist in the general spectral TBA. The properties 
we observed indicate that this function is an important object in general. 

• In this paper we concentrated on the vacuum energy of the flux between the quarks. 
Introducing excitations may lead to an interesting interaction between the conventional 
Bethe Ansatz and the matrix model. 

• Finally, it would be interesting to solve the cusp TBA in other special limits of the 
parameters. For example, in the ladders limit where i9 — > 00 with g e %d l 2 fixed, the 
TBA should reduce to a Schrodinger problem [44]. 

We hope to come back to these points in the future. 

We thank S. Giombi for discussion on this point. 
! We thank N. Drukker for discussion on this point. 
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A Kernels conventions 

For the TBA in section 2 we use the standard definitions 

f [a] (u) = f(u + ia/2) , s(u) = — — r , I m . n = 5 m+hn + <5 m _i. n . (152) 

2 COSh(7TMj 

Our conventions for the kernels are 

r=l 

= \ (K+ - K- + K+ + K-) (153) 

k ba = nt ] + Ct = \ (i^ 6 - 11 - + + + ± K M+2r 

r=l 

where K a , K a and K a are given in (14). The functions 1Z and B are the standard ones 
defined in [8] and are not used in the main text. 



B Derivation of q 

In this appendix we give some more details about the function q defined in section 4.4.2 and 
used in section 4.5.1. 

The function y L+1 Q(y) is defined as the sum of positive powers of y in the Laurent 
expansion of P^y) sinh(27ru), see (79), (84). It is therefore almost proportional to -Pi(y). 
That is, we have 

y L+1 Q(y) = q(y) + P(y)f(y) (154) 
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where q(y) has a finite Laurent expansion and f(y) is some function with a regular expansion. 
One can always shift q by another function proportional to Pjr,. The minimal choice of q is 
such that it has Laurent degree L + l. It is given by 

q(y) = [P(y)smh(2iru)) + -P(y)[sinh(2iru)) + = I te xwnhQiru,) [p{ _ } 

J 2m x z — y z 

It follows from its definition that at the zeros of Pl we have q(yi) = y,f +1 Q(?/j). Moreover, 
for any y it satisfies 

q(l/ y ) = [P(y)sinh(27TM)]_ - (-l) L P(y) [sinh(2vra)]_ (156) 

= (-1) L (P(y) sinh(27TM) - [P(y)wah{2im)] + ) - (-l) L P(y) (sinh(27ra) - [sinh(27TM)] + ) 

= (-i) L+ M</) • 

By this we prove that on the roots of Pl the function y L+1 Q can be replaced by a simpler 
function q(u). 



C Orthogonal Polynomials 

In these appendix we relate the functions Pl to orthogonal polynomials with a simple mea- 
sure. We use this fact to prove that the zeros of Pl{x) are located on the unit circle. This 
was also used in the derivation of a closed expression for the energy in section 4.5.1. 

We defined the coefficients in the Laurent expansion of T(y) as 

Using the explicit form for Pl{v) (95) we can evaluate this integral in terms of the coefficients 

Q 

L+l 

f^Yjd^-L-M- (158) 

»=i 

It follows from (92) f s — for s = 0, . . . , L and Jl+\ = 1- Let us change the variable of 
integration to z = | (y + 1/y) = ^. The integration measure in (157) transforms as 

- + (-D L y s ) ^ = { + ^ dz ' L is even (159) 

y- 1 )y J y \-{2 s z s - 1 + ...)dz , Lis odd 1 ' 

Let us consider separately the case where L is even and odd. 
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Even L. We define T[ +1 (z) = 2irgzP L (x{z)). It is a polynomial of odd degree L + 1 in z 
variable which obeys 



1 

J dz^(z) T°(z) z s = | f ; J = where »(z) = 



sinh(4vr„) ^ 



Angz\Jl — z 2 



At weak coupling g — > the measure reduces to ^=f and becomes a Chebyshev poly- 
nomial of the first kind T n . 



Odd L. Note that for odd L we have P L (±1) = and therefore we can always extract the 
factor y — 1/y — 2\f z 2 — 1. We define U[(z) = 27rg,z J^i^ ■ It is a polynomial of odd degree 
L in z variable which obeys 

i 

[dzfi(z)U^z)z s = l ^ ' S = n where fi(z) = sinh ( 4 ^) ^^z~ 2 , (161) 

-1 

At weak coupling g — > the measure reduces to \/l — -z 2 and become a Chebyshev 
polynomial of the second kind U n . 

The relation to the orthogonal polynomial is quite remarkable because this automatically 
implies that all zeros of Pl(x) can only appear on the unit circle. This follows from the fact 
that all zeros of orthogonal polynomials on the interval [—1,1] must be inside the interval for 
a positive defined measure. Indeed, let's assume that only some of zeroes of an orthogonal 
polynomial p n (z) are real and denote them z n , then the function p n (z) Yl^z — Zj) is either 
always positive or negative on the interval. This is in contradiction with the orthogonality 
which tells that the integral over the interval should vanish. The image of the interval [—1,1] 
under the map \{x + 1/x) — z is the unit circle which doubles the number of roots. 



D Quasi-classical Pl 

Knowledge of the quasimomemtum p(x) and as a result of the resolvent G(x) allows us to 
find the polynomial Pl{x) in the classical limit. In view of importance of this quantity we 
make this calculation here. This calculation mainly consists in the computing of an integral 
of quasimomentum p(x) 

00 00 
g J dx ^1 — — ^ p(x) = —g J dx ^x + — ^ p'(x) — L — g — ^—p(x) (162) 
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where we integrated by part. The reason is that p'(x) is much simpler then p(x) itself and 
does not contain any elliptic functions. As a result the integral can be computed explicitly 



2 / rfx x + - p\x) = y v y . - — y- - AgK + L log A 163 

J \ xj (x 2 -l)y 

X 

L , / 2xV - ^ (2/ 2 + I) 2 + 2y 2 (x 2 - 1) Jy 2 -x 2 J\ - x 2 y 2 

log 1 

2 \ Ax 2 y 2 2x 2 y 

the integral in fact is log-divergent and we introduce a cut-off A. We also denote y = is 
the branch point. At the same time from the relation to the resolvent the integral gives 



<j I <!-r [ 1 - -jp ) /Hi) = log 



' 2L 

n 



.r 



- L log A 



(164) 



and we see that the divergent part cancels. Since log P d = log C\ Y[^ L ^t§ 



and C\ can be 



found in the classical limit to be C\ = exp 



x 2 + 1 



-Llog^iZ^L _ 8 gE 



we finally get 



log P c \x) = -AgE - g p(x) + L log 



2y Ag~K\J y 2 — x 2 a/1 — x 2 y 2 



x 



y 2 -l 



{x 2 -l)y 



(165) 



+ 2^ 



2x A y 2 - x 2 (y 2 + l) z + 2y 2 (x 2 - 1) y^y 2 -x 2 ^l- x 2 y 2 



x 2 (y 2 — l) 2 



+ 2 



x 2 (y 2 - iy 



where again y = e 1 ^. 



E Classical String 

In this appendix we solve the classical string problem at any value of L and the angles 0, 9. 
The solution is given in parametric form in terms of elliptic functions. We then expand the 
general result in various limits. 



E.l Solution for Arbitrary Angles and L 

We are following [15] extending these results to the most general case with arbitrary L, 
and 9. We start from the following ansatz for the embedding coordinates in AdS 3 : 



and S 3 : 



Vi + m = e !K Vl + r 2 (a) , y 3 + iy A = r(a)e 1 ^ 
x x + ix 2 = e !T Vl " P 2 (°) , x 3 + ix 4 = p(a)e if ^ 



(166) 
(167) 
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We use the conformal gauge with energy momentum tensor normalized to 1. In this gauge the 
range of the worldsheet coordinate a G [— s, s] is nontrivial and has to be found dynamically. 
The ends of the strings should go to the boundary in AdS% which implies r(±s/2) = oo. 
Furthermore the angles are defined as <p(±s/2) = ±(n — <p)/2 and /(±s/2) = ±9/2. In order 
for the ansatz to satisfy the equations of motion and the Virasoro constraint, the functions 
0, ip and p, r should satisfy the following equations 



2( J\2 



p\p' 



1-p 2 
1 + r 2 



- P 2 r 

-4 - (7 2 - 1)P 2 + tV - D e 



-£i + (k 2 - l)r 2 + «V = D 



(168) 
(169) 
(170) 
(171) 



From these equations we can write the angles and charges directly, avoiding solving the 
equations themselves. We have 



s = 



s/2 

-s/2 
/2 



da 



dr 



ro 



dp 



2p 



PO 



/S/Z l>OL 
Lp'da = -2 
-s/2 Jr 



dr- 



9 = 



E = 



s/2 
s/2 

-s/2 

4gn 



f'da 



dp- 



PO 



00 , r v / TT72 
dr- 



ro 



finite part 



L = 4#7 / dp 



>y/T 



PO 



(172) 
(173) 
(174) 
(175) 
(176) 



where p (?"o) is a solution to the equation D e = (D^ = 0). The "finite part" means that 
the integral is power-like divergent and one should keep only the finite part. These integrals 
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are standard elliptic integrals and give 



2V2 K f-k 2 d + 1 2 + l 



vV + k 2 + 1 V^ 2 + 7 2 + l 
2L 



2V2 K f ~k 2 +K 2 + l 

K 2 + k 2 + l V K + ^ + 1 



keil + k 2 -^) 
2 7 2 K 



(l + 7 2 + A;, 2 )n 



- 2£| - 7 2 + 1 


k 2 e -l 2 -l 


2k n 


2k n 



k 2 - ry 2 - 1 
2tt^ I 9 7 l 



2ki 



= 



2L 



k^ (1 + k 2 - K 2 ) 

k 2 -K 2 -l 



(1 + k 2 + fcj) n 



fci - 2ft - « a + 1 



2A; 



*2 - k 2 - 1 



2A; 2 



- 2k 2 K 



-2\[2g 



E = -2y/2g 



2k 2 

a/7 2 + k 2 + 1 

7 

/t 2 + k 2 + 1 



E 



E 



-k 2 + 1 2 + l 
k 2 + 7 2 + 1 

-A; 2 + k 2 + 1 
A; 2 + k 2 + 1 



-EC 



- EC 



-A: 2 + 7 2 + l 
A; 2 + 7 2 + 1 

-kl + n 2 + l 

k 2 +K 2 + l 



(177) 



where k e = {Ay 4 - 2 7 2 + 4 7 2 £ 2 + 1 and fc = ^//t 4 - 2/t 2 + 4 7 2 £ 2 + 1. This set of equations 

solves the problem in a parametric form. From the first four equations one should find four 
parameters £g, k and 7 in terms of 0, 6 and L and then the last equation gives the energy. 



E.2 Small angles 

In the small angle limit these equations simplify since then £$ and £g are small as is k — 7. 
Still the solution is given in a parametric form in terms of a new parameter oj defined as 

L = Ag[K(uj 2 ) -E(w 2 )] . (178) 

Let us first consider the case = with 9 small. 

= Case. In this case we see that £$ — and therefore k^ = \ / k, 2 — 1. The equations 
for s and E now reduce to 

s = -K(1/k 2 ) (179) 

K 

E = 4g [K{1/k 2 ) - E(1/k 2 )] . (180) 
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Furthermore, if we expand the reminding equations to the leading order in £g we get 



7 



-K(l/ 7 2 ) - t 2 



( 7 3 + 7 )E(l/ 7 2 ) 7 K(l/ 7 2 ) n 
( 7 2 -l) 2 1"7 2 

2K(l/ 7 2 ) 4 7 2 E(l/ 7 2 ) 



L = 4^[K(l/ 7 2 )-E(l/ 7 2 )] + 
7 Z — i 



1 - 7 2 



+ 0$) (181) 
+ 0(4) (182) 



(7 2 " 1) 



(183) 



We want to express the result in terms of u defined in (178). From (182) we find 7 in terms 
of u 



7 = - + 

UJ 



1 K(u 2 ) 



lo(lo 2 - 1) 2wE(w 2 ) 



Then from (179) and (181) we find « 



1 



K 



Substituting this into (180) and using (183) we get 

„ T „ 2 u; 2 E(u; 2 ) n2 1 - to 2 
E-L = 2g£j - 1 / = fffl 2 



(184) 



(185) 



1 -w 2 



2E(w 2 ) 



(186) 



This is the result used in the main text. 
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General case. The procedure explained in the previous section can be applied to get 
expansion in both <fi and 9. We quote the result only 

/ K 2 (a; 2 - if K(to 2 -9) (to 2 -if (a; 6 -7uj 2 + 6) \ 
^ 32E 5 o; 2 + 96E 4 o; 2 + 96E 3 o; 2 J 



E-L 1-to 2 

+ e 2 



g((j> 2 - # 2 ) 2E 



K 2 (o; 2 -l) 3 K(5o; 2 + 3)(a; 2 -l) 2 (5a; 4 - 6a; 2 + 1 A 
^ ^ 32E 5 o; 2 96E 4 o; 2 96E 3 J { ' 

4 / 7K 4 (a; 2 - l) 5 _ K 3 (9a; 2 - 37) (a; 2 - if K 2 (17a; 4 - 138a; 2 + 225) (a; 2 - if 
+ I 768E 9 o; 4 768E 8 o; 4 + 2304E 7 o; 4 

K (2a; 6 + 281a; 4 - 1110a; 2 + 1035) (to 2 - if 
+ 11520E 6 o; 4 

(to 2 - 1) (to 8 - 16a; 6 + 241a; 4 - 570a; 2 + 360) \ 
+ 11520E 5 w 4 J 

2 2 / 7K 4 (to 2 - if _ K 3 (3a; 2 - 17) (to 2 - if _ K 2 (2a; 4 + 27a; 2 - 81) (a; 2 - l) 3 
+ ^ ^ 768E 9 a; 4 384E 8 a; 4 1152E 7 a; 4 

K (14a; 6 - 28a; 4 + 135a; 2 - 225) (a; 2 - l) 2 
5760E 6 a; 4 

(a; 2 - 1) (-14a; 8 + 14a; 6 + 31a; 4 - 60a; 2 + 45) \ 
+ 11520E 5 a; 4 J 

4 / 7K 4 (a; 2 - l) 5 K 3 (9a; 2 + 19) (a; 2 - l) 4 K 2 (29a; 4 + 30a; 2 + 45) (a; 2 - l) 3 
+< P [ 768E 9 a; 4 + 768E 8 a; 4 + 2304E 7 a; 4 

K (122a; 6 + 11a; 4 + 30a; 2 + 45) (a; 2 - l) 2 (a; 2 - 1) (61a; 4 - 46a; 2 + 1) \ 
11520E 6 a; 4 + 11520E 5 J 



+ 



Here IK = IK(a; 2 ) and similar for E. It is easy to get more orders in and 9, however, the 
expressions become very bulky. We found them explicitly up to the order 10 and can provide 
them by request. 
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Small L. Further expanding at small L, we find 

E-L f 1 O 2 - 5<f) 2 -6 A + U<f) 2 6 2 - 37(f) 4 . 

+ ^ (188) 



(f) 2 - 6 2 » \ 7T 8tt 3 64tt 5 

3 _ 3 (6 2 - 17(f) 2 ) + 3 (5# 4 - 140 2 # 2 + 730 4 



+ 



4tt 2 64vr 4 256tt 6 

L 2 ( 9 3 (2fl 2 - 90 2 ) 3 (13fl 4 - 1380 2 fl 2 + 6530 4 



g V 64tt 3 128vr 5 8192tt 7 



L 3 f 5 -133# 2 - 3790 2 -41# 4 - 20 2 # 2 - 3410 



+ 



g 2 V 256vr 4 8192tt 6 4096tt 8 

L 4 / 45 105(11# 2 - 70 2 ) 3 (2OO10 4 + 36980 2 fl 2 - 1390 4 



g 3 V16384tt 5 131072tt 7 1048576tt 9 

The first line reproduces the known result of [45]. For 9 ~ we get 

E-L ( g 3L 9L 2 5L 3 45L 4 \ , , 

/ + 1CT ,..uJ ( 18 9) 



9 2 \ n 4tt 2 64#7r 256# 2 tt 16384c/% 5 

o 3L 21L 2 L 3 105L 4 

+ 4 1 - A + 



2tt 3 4tt 4 128#tt 16# 2 tt 32768#% 
3g 3L 99L 2 3L 3 2085L 4 



3vr 5 4tt 6 512#7r 32# 2 tt 131072# 3 7r 

5g 3L 225L 2 L 3 7905L 4 
+ 0" I + 



16tt 7 4tt 8 1024c/tt 8# 2 tt 262144c/ 3 tt 
35o 3L 1995L 2 5L 3 97425L 4 

+ 0° I -TTrA + 



128vr 9 4tt 10 8192^7r n 32# 2 tt 20971520% 
Large L. For L/g to be large, u should be exponentially close to 1 

U = 1 + M + f (32- ^) + ^( 4 "^ ~ + ^) (190) 

+ ^_5^ + 38^_134« N 
V 3# 3 g 2 g J 

— — —2 

where 5 = e 2 a . One can plug this expansion of uj into the expansions (187) to get the 
large L expansion. We do not use this limit in the paper so we skip the details. 
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F Strong Coupling Expansion for Fixed L 



Here we collect the strong coupling expansion of our result for various L's. 

L B L + 0{l/g^ 



n £ _ _(_ a IT 

U vr 8vr 2 ^ 128tt 3 c/ ^ 

1 £ 9__i_ 39 

1 - 8tt 2 "t" 128tt 3 c/ "T" 



a £ 27 |_ 363 , 

vr 8tt 2 ~ r 1287r 3 o 

c; £ _ _33_ i 543 , 

° a T 2 T 1o0 „3„ I 



3 


+ 
+ 


63 


5127r 4 g 2 
39 


32768ti-V 
279 


5127r 4 o 2 
165 


32768tt 5 o 3 
1449 


512tt 4 o 2 
441 




32768ttV 
9441 


5127r 4 o 2 
927 




327687T V 
30177 


512tt 4 o 2 
1683 




32768tt 5 g 3 
72297 


512tt 4 9 2 




327687T 5 g 3 



o 9 15 I 111 I 100 I44jj Mrm 

Z vr 8vr 2 128tt 3 c/ ^ 512tt 4 9 2 327687r 5 .g 3 

Q 9_ _2J_ i 219 , 

d - 8vr 2 ^ 128vr 3 c, ^ 



128tT 3 £ 

This table can be summarized by the following single equation 

g 6L + 3 3(6L 2 + 6L + 1) 10L 3 + 15L 2 + 11L + 3 



£ 7 



vr 8tt 2 128vr 3 # 512ttV 

9 (10L 4 + 20L 3 - 22L 2 - 32L - 7) 



32768ttV 



(9(1//) . (192) 



We see that knowing just a few first L's we can in fact extrapolate to large L using the 
polynomial property of the coefficients. In particular we can make a comparison with the 
classical limit by re-expanding for large L ~ g a la [46]. Denoting £ = L/g we get 

/l 3£ 9£ 2 5£ 3 45£ 4 m ,~\ , . 

o / 3 9£ 15£ 2 45£ 3 , 4 v 

+ 9 V 8^ + 647T 3 + 5127T 4 ~ 81927T 5 + ' ' 

1/3 11£ 99£ 2 , 3 v 

+ # Vl28^3 + 5127 4 + 16384tt 5 + ^ > 



1/3 9£ 



# 2 V5127T 4 10247T 



+o ( £2 )) + M» + °( £i )) + °( 1 ^ 



we see that the first line by dimension is the classical contribution. The second line is the 
one-loop part of the result and so on. Indeed, this procedure gives the correct classical energy 
in the first line of (189)! This is another more direct test of our results which does not involve 
the algebraic curve consideration. Note that the first term in the second line in the above 
equation perfectly reproduces the one loop result of [45]. 
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